Ground state properties of a confined simple atom by Ceo 

fullerene 



M. Neek-Amal, 1 G. Tayebirad, 2 and R. Asgari * 3 

Department of Nano-Science, Institute for Studies 
in Theoretical Physics and Mathematics (IPM), 
P.O. Box 19395- 5531, Tehran, Iran. 
2 Department of Physics, Iran University of Science and Technology, Tehran, Iran. 
3 Institute for Studies in Theoretical Physics and Mathematics, Tehran 19395-5531, Iran. 

Abstract 

We numerically study the ground state properties of endohedrally confined hydrogen (H) or 
helium (He) atom by a molecule of Cgo- Our study is based on Diffusion Monte Carlo method. We 
calculate the effects of centered and small off-centered H- or He-atom on the ground state properties 
of the systems and describe the variation of ground state energies due to the Cgo parameters and 
the confined atomic nuclei positions. Finally we calculate the electron distributions in x — z plane 
in a wide range of Cgo parameters. 
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1. INTRODUCTION 



The Properties of Carbon allotropes with closed cage structures have been an active sub- 
ject of research since its first experimental discovery back in 1985 Accurate calculation of 
the structural properties of single shell- or tubes fullerenes have been performed by Haser et 
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al. |2J using the Hartree-Fock methods, by Dunlap et al. [3j and Heggie et al. |4| using density 
functional theory. The later group considered fullerenes up to a large number of atoms. The 
unique spherically shaped cage of a Ceo fullerene contains empty space that is large enough 
to incorporate atoms of any kind of endohedral composite. The molecule is not very reactive 
since the carbon atoms are all completely saturated by taking part in an extended it— bonded 
electron system similar to graphite. The interior of Cm is expected to have a good capacity 
for confining small molecules, ions and atoms j^, [f], [z, 8, 9|. Physical properties of a con- 
fined simple atom by Cgo fullerene is a subject of significant interest in recent years, mostly 
because they exhibit properties that can lead to important applications in nanostructure 
science and technology 10j. Furthermore, many researchers have been working on modeling 
the storage of various types of gases and metal atom in nanotubes and fullerenes [ll, \\\ . 
This type of research is also fundamental in the application of nano-technology to gas stor- 
age devices, used in such systems as fuel cells, construction of molecular sieves and filtration 
membranes. 

In addition to the real atomic interaction approach, one can modeled the interaction 
between adsorbed atoms to Ceo molecule by an effective attractive potential which can hold 
electrons of atom near Cm wall. This effective potential, in effect, replaces the C m w alls and 
has been successful in explaining many experimental and computational results [kI \l\ . It 
is physically expected to observe significant difference between the ground state properties 
of the confined atoms by Ceo fullerene, and their free states [ll]]. Obviously, in this case the 
strength of the attractive potential is not so strong as to immobilize electrons. Recently, 
some authors 

HQ 

used a similar model of Cgo molecule constructed by Xu et al. [14J to 
calculate the properties of the endohedrally confined atoms. 

Researches using modeling studies in this field have employed various computational 
methods based on prescribed inter-atomic potential energy functions. The Schrodinger 
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equation provides the accepted description for microscopic phenomena at non-relativistic 
energies. This equation can be solved analytically only in a few highly idealized models. 
Calculation of the ground state properties of a system with several atoms/molecules or with 
an arbitrary boundary conditions is a scientific efforts. A theory based on Quantum Monte 
Carlo (QMC) methods allows us to calculate the ground state properties more accurately. 
For the many body system, QMC has been very useful in providing exact results. 

The main purpose of the present work is to calculate the effect of a centered or small 
off-centered Hydrogen (H) or Helium (He) atom on the ground state energy and describe 
also the variation of the ground state energies on the strength of attractive model potential 
within diffusion Monte Carlo (DMC) method. Although the simple spherical shell model for 
Cqo molecule does not completely describe the real endohedrally off-centered atoms, but we 
use it for considering the off-centered effects as a first approximation. 

The contents of the paper are briefly as follows. In section II, we review the basic 
formalism, including our model to the system and in the next section we briefly discuss the 
DMC method with the appropriate steps ,as well as the employed trial wave function. Our 
numerical results for the atoms confined by Ceo fullerene are presented in Sec. IV. Finally, 
in Sec. V we summarize our main conclusion and discuss our further works. 

2. THEORY 

Our approach is based on modeling an endohedral medium by a short-range spherical 
shell with an attractive potential Uq in the range a <r < ka, where a and ka are the inner 
and outer radii of the shell with A = a(k — 1) as the thickness of the shell. Fig. 1 depicts 
all the information on our selected model in the case of a H-atom confined endohedrally 
within a C 60 fullerene. In the case of He-atom we have a similar model and a figure with 
two electrons can be fancied. 

When an electron associated with an atom enters the shell medium, it will be affected by 
further attractive potential, Uq, in addition to the Columbic potentials due to other electrons 
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and nuclei. Therefore the potential energy is written as 

(E^pp^-^ otherwise, 
where x is the position of nucleus, f\ represents the position of ith electron which is 
measured from geometrical center of the shell and Z is the atomic number. Summation runs 
over the number of electrons. The atomic nucleus is always assumed to be fixed at position 
x. Moreover, we calculated the effect of small off-centered nuclear position on ground state 
properties by moving the position of nucleus along a direction and assume it fixed in each 
simulation. 

I — I 

Xu et a/.[14| found a semi-empirical function for the potential of Ceo fullerene and cal- 
culated the parameters a and b by fitting them to experimental data. Specifically, they 
found a « 5.75 au which is approximately equal to the radius of Ceo an d A « 1.89 au. 
The value of the shell thickness, A, is assumed to be different by different authors [13]. We 
fixed a and A by the above values in most parts of our calculations but in order to take 
the variation of A into account, we also calculate the ground state energy as a function of 
thickness by considering that for every simulation these parameters are fixed. Although the 
effective value of the attractive potential Uo is used in Ref. from 0.3015 to 0.4228 au, 
we examined our results in a wide range of values for Uq. This is because of the number of 
carbon atoms which can be changed in different systems. 



19, 



20 



211 ] to calculate the ground state 



We have used the well-known DMC method 
energies, E of simple endohedrally confined atoms such as H- or He-atom. The number 
of electrons for such a system are less than three and thus the ground state wave function 



has no node, consequently the Pauli exclusion princip 



is necessary to use the fixed-node QMC calculation 



e does not play a role. Note that it 



20j or Exact Cancelation method for a 



system with large number of fermions |22|. In what follows, we briefly discuss the DMC 
method used in our numerical calculations. 

Through an analytic continuation of time t to imaginary values — it, the Schrodinger 
equation of the system is given by 

m £ T) = DV 2 MR, r) + {E T - V(R))i>(R, t) , (2) 

4 



where the constant D = H/2m, m is electron mass and Et is an offset energy which 
can be adjusted at the beginning from contributing an average potential for all the initial 
positions of particles. In the limit r — > oo, ?/>(R, r) tends to either a non-zero value or an 
infinite value provided by Et- 

To avoid the effect of singularity in the potential, it would be essential to use a guiding 
trial wave function. More specifically, at the positions where ^(R) diverges, the creation or 
annihilation rate, Et—V(R) becomes tremendously large and makes the numerical algorithm 
unstable. In this case, we used a time independent guide function ^t(R) by introducing a 
function p(R, r) giving by p(R, r) = ^ T (R)^(R, r). By substituting p(R, r) in Eq. ([2]) we 
have 

= DV R [V R -F(R)]KR,r) 

- (E l (R)-Et) P (R,t) , (3) 

where -F(R) = 2V?/ ; t(R)/V ; t(R) is the so-called Fokker- Planck force and corresponds to 
a force function which drifts or walks away from regions where |?/>t(R)| 2 becomes small. 
Meanwhile, the local energy function is El(R)= — -DVr?/>t(R)/Vt(R) + ^(R)- It is impor- 
tant to choose the guide function such that -Ex(R) does not have any singularity. We used 
the Green function method to solve Eq. (J3]). We follow the standard way to decompose the 
Green function into two parts. The usual diffusive part is given by 22] 



/ (R/— R— P-F(R)At) 2 n 

G D (R, R, Ar) = — ^ L + o(Ar)2 , (4 ) 
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(4tt D At 

and the branching term which is given by 

G B (R, R, Ar) = exp[-Ar( ^ (R) + ^ (R/) - E T )) . (5) 

In order to reduce numerical errors to order (Ar) 2 , we have used the Metropolis procedure. 
After performing the diffusion and branching processes for all the particles, the final MC 
step lies in updating the value of Et on an ensemble average. 
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where Nq refers to the initial number of walkers and N denotes number of walkers updated. 
This sort of adapting is essential to avoid large amount of fluctuations in the number 
of walkers. The ground state of the system is obtained by averaging Et over many MC 
steps. Prior to presenting our results, it would be illustrative to discuss the numerical errors 
in the DMC method. Principally, there are two types of errors which limits the accuracy 
of most DMC calculations: (a) Statistical or sampling errors associated with the limited 
number of independent sample energies used in determining the ground state energy, (b) The 
systematic errors associated with finite time-step At, round-off in computing, imperfectness 
of random number generators etc. The total errors in our numerical calculation is about 
for more or less all cases. 

3. NUMERICAL RESULTS 

We started by focusing on the centered H-atom. Our choice of radial part of ground state 
wave function is the linear composition of an unperturbed H-atom wave function and outer 
well wave function. For having the unique function with respect to Uq, we calculated the 
radial ground state wave function by creating the same results as in Ref. jl5j for several 
values of Uq and then we found the best fit with some constants depending to Uq value. Our 
trial wave function for every Uq value is 



where tq is approximated by the shell's mean point and is about a + A/2. 

When Uq is close to zero, we have a free electron ground state properties and when Uq 
is large enough, we have square cosine function for a wave function. By using the cusp 
condition which leads to ip' H {r) + ipH{ r ) — 0, we can determine parameter B in terms of 
other parameters as 



This condition guarantees to cancel singularity at small r— region. From the fitting process 



^H{f) = Aexp(-Br) + C exp(—E(r — tq) 2 ) , 



(7) 



B 



qMO) + 2CEr exp(-Erl 
A 



(8) 
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mentioned above, we obtained 



A = 3.06exp(2.99 U ) 

_ 0.146 Uo 2 -0.03 Up (Q) 
£/ 2 -0.58 C/o+0.33 

E 



1.34 t/ 2 -0.17 C/p 
l/ 2 -2.02 f/o+0.81 



For the He case in our system, we build a suitable trial function which is a product of 
two simple H-atom trial wave function connected by a Pade-Jastrow function 



where r 12 



i'He(r 1 ,r 2 ) = if)H{n)i(>H(r2)ip(ri2), (10) 
f\ — r^l, 9?(r 12 ) = exp(r 12 /2(l + qt 12 )) and a = 0.2 is a well-known value for 



He system 23( . The parameter B is fixed by using the cusp condition for the He case 



2^(0) + 2CEroexp(-Er 2 ) 
5 = ^ • (11) 

Note that the pre-factor 2 in the first term of the nominator comes from the fact that there 
are two electrons in the system with respect to the H system. 

To Consider small values of off-center ( 8 < a) for both H and He cases, we have assumed 
that the wave function of the system is again given by Eq. (JTj) and Eq. fTTUT) and one only 
needs to shift the r parameter by a small amount so that r — * f+5 where 5 is small. 
Furthermore we restrict our calculations to a range such that S becomes less than 0.5a along 
the z-axis. 

In our simulations we fixed At m 0.05. ( See the appendix for more details). Throughout 
the paper, we have implemented the following atomic units. Here and hereafter we use the 
Bohr radius, h 2 /me 2 = 0.529 A as the length unit, the value of h 3 /me 4 = 2.419 x 10~ 2 fs as 
time unit and the Hartree, me 4 /h 2 = 27.25 eV as the energy unit. The main results of our 
work are shown in Figures 2-7. 



3.1. H-atom endohedrally confined in Cgo molecule 

First of all we calculated the ground state energy of H-atom endohedrally confined within 
Ceo molecule as a function of Uq, the confining well strength for different small nuclei posi- 
tions. The nucleus vector position is considered as x = (0, 0, Da) where D assumes different 
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values of 0, 0.25 and 0.5. For off-centered cases, H-atom obviously does not have any spheri- 
cal symmetry and as expected from an analytical approach, it becomes either more difficult 
or undoable for such configurations. It is apparently obvious from Fig. 2a that for the small 
[/p-values, the ground state energy remains constant at the usual well-known value of -0.5 
for ground state energy (-13.6 eV). It means that the electron is not influenced by the outer 
confining well. For the case of D = 0, our results are in excellent agreement with recent ex- 
act results which were obtained by solving the Schrodinger equation directly From this 
figure, when the H-atom moves towards to the confining well region, the energy maintains 
it's -0.5 value. It then starts to decrease, when it reaches the threshold point at Uq ~ —0.7. 
The value of Uq at the threshold point depends on the assumed values of D. As physically 
expected by moving the nucleus of the H-atom far from the shell, the threshold value of Uq 
should move to the right because the atom approaches the free atom case. The effect of 
thickness for the ground state energy is shown in Fig. lb with a = 5.75. The electron of the 
confined H-atom becomes more immobile when the shell thickness increases. Consequently 
the electron becomes localized by increasing the shell thickness or the value of the attractive 
potential. For sake of clarification, the numerical data of the ground state energies are given 
in the Tables I and II. 

The probability of finding the electron at a distance r from its nucleus is proportional to 
r 2 p(r) where p(r) = (p{r)i(jH{f) and 0(r) is the electron radial wave function assuming that 
the nucleus of the H-atom remains at the origin. Fig. 3 depicts the r<fi{r) as a function of r 
for different f/o-values. The electron does not feel the effect of the well for small [7o-values, 
however it is trapped in the shell region for large Uo-va\xies. It is important to note that the 
r0(r) does not have any tail for very large attractive potential. 

When the H-atom is off-centered, we calculated two-dimensional electron distributions in 
the x — z plane. Fig. 4 shows the effects of centered (D = 0) and off-centered (D = 0.25) 
H-atom on the electron distributions at threshold value of Uq. In the symmetrical case, 
electron distribution is uniformly distributed. Our results show that electron distribution 
has a peak around the nucleus in all directions for a weak attractive potential. However, 
the electron is localized in the vicinity of the shell region for a strong attractive potential. 
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In this case, the atom becomes ionized depending on the large value of Uq and D. As it is 
physically understandable, the electron distributions along x— or y-axis holds its symmetry. 
However the electron distribution along z-axis has no symmetry. 

3.2. He-atom endohedrally confined in Cqq molecule 

To the best of our knowledge, there is no analytical solution for He-atom confined within 
Ceo molecules. Recently, some authors have intensively studied the two-electron photoion- 
ization cross section of He atom confined in Cgo-pj] They calculated the regular and irregular 
solutions of the Schrodinger equation for the system. Here we present more accurate ground 
state energy as well as the associate wave function for the system which will be used to 
calculate other physical quantities. 

We considered a He-atom bounded by fullerene and calculated the ground state proper- 
ties. The He-atom with its two electrons and nucleus centered at x has an anti-symmetric 
excited state wave function. The method described in this paper allowed us to calculate the 
ground state energy as functions of attractive potential and the well thickness. 

In Fig. 5 we show the ground state energies as a function of attractive energy, Uq for 
different D and thickness values. Similar to the H-atom results, electrons become more 
immobile with increasing the shell thickness. Note that in the limit of zero Uq, the ground 
state energy gets the well-known value of the He ground state energy which is —2.903 
(— 78.96eV). The ground state of energy decreases after the threshold Uq ~ —1.9 value for 
the centered case and the threshold value is difference from the H-atom problem. Numerical 
values of ground state energies are given in the Table IV. 

Finally, in Fig. 6 we show electron distributions in the x — z plane when nucleus remains 
at origin for a high value of Uq. Clearly, two electrons prefer to stay far from each other. 
Moreover by increasing the value of Uq, the probability of finding each electron near the 
shell increases. 
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4. CONCLUDING REMARKS 



In this work we have numerically calculated the ground state energies and electron distri- 
butions of a simple atom confined by Ceo fullerene within the diffusion Monte Carlo method. 
We modeled the Ceo fullerene with attractive confining well and obtained the physical quan- 
tities of one or two electrons of the atom endohedrally. Our results for a hydrogen atom 
located at the origin is in excellent agreement with exact calculations. Within the DMC 
method we obtained the ground state energy for off-centered H- and He-atoms for different 
confining well potentials. We described the electrons distributions for endohedral case as 
a function of nucleus positions and the confining well potential as well. Moreover, elec- 
tron distributions are more influenced by changing the nucleus position. We believe that 
for an accurate quantitative calculations of off-centered effects, the real physical model for 
Ceo molecule should be deformed in the presence of endohedral atoms and deviation from 
spherical symmetry should be considered. 
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APPENDIX A: DETAILS ON EXPLICIT CALCULATION OF THE TIME STEP 

The time step can be selected according to the criterion that the systematic error asso- 
ciated with the use of a finite Ar is less than the statistical error. To achieve the desired 
variance in energy, the value of N and run time are chosen sufficiently large. In these 
simulations we chose N and the run time in the range of 4000-6000. In Fig 7. we show the 
energy and its standard deviation as a function of Ar for free He-atom. As it is clear from 
the figure, the range [0.05 — 0.1] should be the acceptable range of the time step, Ar. 
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FIG. 1: (Color online) The schematic representation of H-atom inside an attractive shell for mod- 
eling the H-atom endohedrally confined inside Ceo molecules. 
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FIG. 2: (Color online) Top: Ground state energy for H-atom located at position x = (0, 0, Da) as 
a function of Uq in units of Hartree. Bottom: Ground state energy of H-atom located at origin as 
a function of Uq in units of Hartree for different thickness values. Here a is 5.75. 
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FIG. 3: (Color online) Radial electron wave function multiply to position, r<fi(r) as a function of r 
for H-atom located at D = for different Uq values. Dashed vertical lines refer to the shell region. 
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FIG. 4: (Color online) Electron distributions of H-atom in x — z plane for Uq = —0.7 where D 
(top) and D = 0.25 (bottom). Legends denote the number of normalized electron density. 
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FIG. 5: (Color online) Top: Ground state energy of He-atom as a function of Uq in units of Hartree. 
The nuclear is located at x = (0, 0, Da) for D = 0, 0.25 and 0.5. Bottom: Ground state energy of 
He-atom as a function of Uq in units of Hartree for different thickness values. Here a is 5.75. 
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FIG. 6: (Color online) Electron distributions in x — z plane for Uq = —4.0 for two electrons of 
He-atom confined at D = 0.0. Legends denote the number of normalized electron density. 
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TABLE I: Ground state energy for H-atom for different D values. The statistical error bar for 
GDMC is about ±10" 4 . 
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TABLE II: Ground state energy for H-atom for different thickness values. 
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D=0.0 


D=0.25 


D=0.5 


0.00 


-2.9031 


-2.9036 


-2.9038 


0.20 


-2.9033 


-2.9045 


-2.9103 


0.40 


-2.9040 


-2.9057 


-2.9116 


0.60 


-2.9044 


-2.9063 


-2.9138 


0.80 


-2.9055 


-2.9077 


-2.9144 


1.00 


-2.9061 


-2.9088 


-2.9152 


1.20 


-2.9064 


-2.9145 


-2.9169 


1.40 


-2.9073 


-2.9188 


-2.9176 


1.60 


-2.9081 


-2.9231 


-2.9276 


1.80 


-2.9087 


-2.9424 


-2.9445 


2.00 


-3.0523 


-3.0560 


-3.0653 


2.20 


-3.4101 


-3.4148 


-3.6131 


2.40 


-3.7751 


-3.7756 


-3.9799 


2.60 


-4.1364 


-4.1441 


-4.3421 


2.80 


-4.5043 


-4.5053 


-4.7079 


3.00 


-4.8758 


-4.8705 


-5.0756 


3.20 


-5.2496 


-5.2502 


-5.4505 


3.40 


-5.6228 


-5.6228 


-5.8193 


3.60 


-5.9986 


-5.9976 


-6.1959 


3.80 


-6.3742 


-6.3767 


-6.5741 


4.00 


-6.7538 


-6.7494 


-6.9580 



TABLE III: Ground state energy for He-atom for different D values. 
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-U 


A = 1.89 


A = 3.745 


A = 5.6 


0.00 


-2.9031 


-2.9032 


-2.9032 


0.20 


-2.9252 


-2.9258 


-2.9263 


0.40 


-2.9260 


-2.9261 


-2.9265 


0.60 


-2.9333 


-2.9363 


-2.9372 


0.80 


-2.9344 


-2.9372 


-2.9394 


1.00 


-2.9388 


-2.9396 


-2.9429 


1.20 


-2.9435 


-2.9455 


-2.9468 


1.40 


-2.9453 


-2.9469 


-2.9479 


1.60 


-2.9559 


-2.9614 


-2.9738 


1.80 


-2.9657 


-3.2565 


-3.2902 


2.00 


-3.0823 


-3.6351 


-3.6563 


2.20 


-3.4101 


-4.0522 


-4.0453 


2.40 


-3.7751 


-4.4060 


-4.3871 


2.60 


-4.1364 


-4.8250 


-4.8418 


2.80 


-4.5043 


-5.1915 


-5.1753 


3.00 


-4.8758 


-5.5377 


-5.5693 


3.20 


-5.2496 


-5.8922 


-5.9155 


3.40 


-5.6228 


-6.3206 


-6.2839 


3.60 


-5.9986 


-6.6892 


-6.6714 


3.80 


-6.3742 


-7.0863 


-7.0453 


4.00 


-6.7538 


-7.4316 


-7.4289 



TABLE IV: Ground state energy for He-atom for different thickness values. 
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